Time evolution of the Partridge-Barton 

Model 

Roberto N. Onody*and Nazareno G. F. de Medeiros ' 

Departamento de Ftsica e Informdtica 
Instituto de Ftsica de Sao Carlos 
Universidade de Sao Paulo - Caixa Postal 369 
13560-970 - Sao Carlos, Sdo Paulo, Brasil. 



Abstract 

The time evolution of the Partridge-Barton model in the presence of 
the pleiotropic constraint and deleterious somatic mutations is exactly 
solved for arbitrary fecundity in the context of a matricial formalism. 
Analytical expressions for the time dependence of the mean survival pro- 
babilities are derived. Using the fact that the asymptotic behavior for 
large time t is controlled by the largest matrix eigenvalue, we obtain the 
steady state values for the mean survival probabilities and the Malthusian 
growth exponent. The mean age of the population exhibits a t -1 power 
law decayment. Some Monte Carlo simulations were also performed and 
they corroborated our theoretical results. 
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1 Introduction 



Early in life we perceive that everything around us, inanimate objects, ani- 
mals and human beings undergo a variety of changes that accompany the pas- 
sage of time. Everything suffers a progressive deterioration with time. This 
phenomenon is called ageing or senescence and it is characterized by a decline 
in the physical capabilities of the individuals. Several theories (see |jj and ref- 
erences therein) have been suggested to explain why there is senescence, when 
it occurs and what are the biological processes responsible for it. Usually, these 
theories are divided into three classes: biochemical, evolutionary and telom- 
eric. The first invokes damages on DNA, cells, tissues and organs and connect 
senescence with imperfections of the biochemical processes. One kind of this 
biochemical imperfection is the presence of free radicals which can cause death 
of the cells or may even lead to cancer [|| . The evolutionary theory || || , on the 
other hand, explains the senescence as a competitive result of the reproductive 
rate, mutation, heredity and natural selection. In the telomere hypothesis ||, 
senescence depends on the cumulative number of cell divisions. The replication 
of a normal cell is followed by a tclomcric shortening. This acts as a counting 
mechanism which controls the number of divisions. 

Evolutionary theories of ageing are hypothetico-deductive in character, not 
inductive. They do not contain any specific genetic parameter, but only phys- 
iological factors and constraints imposed by the environment. There are two 
kinds: the optimality theory and the mutational theory. In the optimality the- 
ory H , senescence is a result of searching an optimal life history where survival 
late in life is sacrificed for the sake of early reproduction. For the mutational 
theory (J, , on the other hand, ageing is a process which comes from a balance 
between Darwinian selection and accumulation of mutations. The natural selec- 
tion efficiency to remove harmful alleles in a population depends on when in the 
lifespan they come to express. Alleles responsible for lethal diseases that express 
late in life, escape from the natural selection and accumulate in the population, 
provoking senescence. Nevertheless, if the natural selection is too strong then 
deleterious mutations might not accumulate in the population and the eternal 
youth could be reached. An evolutionary model with such characteristics was 
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recently studied and solved by Onody and de Medeiros M. 

A simple evolutionary model of ageing is the Partridge-Barton model||. 
It was introduced to illustrate the optimality theories of ageing. Its principal 
feature is the inclusion of the antagonistic balance mechanism JlO| . This mecha- 
nism arises out from processes which enhance the lifespan early in life, but have 
deleterious effects latter. 

In this work, we find an exact solution for the whole dynamics of the 
Partridge-Barton model. When only deleterious somatic mutations and pleiotropy 
are present the time evolution of the model can be formulated in a matricial 
form. Explicit analytic expressions can be written for the mean survival proba- 
bilities and the growth rate. For large time t, the system behavior is dominated 
by the matrix largest eigenvalue. The existent integrals can be solved by the 
saddle point approximation, allowing us to determine precisely the steady state 
values of the survival probabilities. A time expansion for the population's mean 
age shows that it converges to a constant value according to at -1 power law, a 
result which was first obtained by Ray All the results were confirmed by 
some Monte Carlo simulations that we performed. 

2 The Partridge-Barton Model 

In the Partridge-Barton model there are only three ages. The population 
consists of babies (age = 0), juveniles (age = 1) and adults (age = 2). The sur- 
vival probabilities from infancy to juvenile is J\ and from juvenile to adulthood 
is J2. Reproduction is permitted only to juveniles and adults, with rates mi 
and m2, respectively. Babies don't have offsprings and adults are eliminated 
from the population after reproduction. 

The population grows at a steady rate r. The Malthusian growth exponent 
r is related to the other parameters of the model through a discrete version of 
the Euler-Lotka equation ]T^j 

777.1 J\e~ r + m2J\Jie~' lr = 1. (1) 

The antagonistic pleiotropy fllPf arises when the same gene is responsible 
for multiple effects. For example, genes enhancing early survival by promotion 
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of bone hardening might reduce later survival by promoting arterial hardening. 
Partridge and Barton implemented the basic ideas of the antagonistic pleiotropy 
by adopting the constraint, J\ + Jf = 1, between the survival probabilities Ji 
and J2. The parameter x is a real positive number whose value depends on the 
kind of population we are dealing with. The pleiotropic condition states that it 
is impossible to sustain simultaneously both high juvenile and adult survivals. 
For the particular case in which mi = ra-i — 1 and x = 4, Partridge and Barton 
found Ji = 0.935 and J2 = 0.505 as the values which maximize the growth rate 
r. 

Also the action of deleterious or helpful mutations can be added to the 
model. Using Monte Carlo simulations, Stauffer studied the case in which 
the pleitropic constraint J% + J% = 1 is accompanied by random somatic muta- 
tions. His results clearly show that the survival probabilities J\ and J2 move 
rapidly to stationary values with J\ > Ji. This fact means that the model ex- 
hibits senescence, in the sense that the adult survival is lower than the juvenile. 
In the absence of mutations, J\ and J2 tend towards 0.935 and 0.505 in accord 
with the Partridge-Barton conclusions. However, it is not clear how the system 
drives itself towards these optimal values. 

3 Analytical Solution 

In this section we obtain the exact time solution of the Partridge-Barton 
model in the presence of pleiotropy and somatic mutations. 

Let Ni(Ji, t) be the number of individuals at age i (i = 0, 1, 2) with survival 
probability between Ji and Ji + dJi at time t. We choose, as initial condition, 
a population with the profile 

Nt(Ji,0) = N Q S ifi , (2) 

that is, in t = there are only iVo babies with the survival probabilities Jo 
uniformly distributed in the interval [0, 1]. 

At time t, all babies are equally submitted to somatic and deleterious mu- 
tations with strength a (a < 1). Their survival probabilities Jo are changed 
to J\ — olJq. Subsequently, all these babies pass through natural selection in a 
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such way that, on average, the number of juveniles with survival probability J\ 
at the instant t + 1 is given by 

N 1 (J u t+l) = J 1 N Q (Jo,t). (3) 

Since the mutation is restricted to be somatic, each one of the Aq(Ji,i + 1) 
juveniles will give birth to exactly mi offspring with survival probability J . 

Now, the probability with which a juvenile will reach adulthood must take 
into account the antagonistic pleiotropy and the somatic deleterious mutations. 
As pleiotropy is not affected by the somatic mutations, a juvenile with survival 
probability Ji (formerly, a baby with survival probability Jo) will change its 
survival probability to (f — Jo) 1//x , where x is a real positive number and a 
measurement of the pleiotropic constraint. Under the action of a deleterious 
somatic mutation, described by a parameter (0 < 1, fixed), the new survival 
probability can be written as J 2 = 0(1 — Jo) 1 ^- Submitting all juveniles to nat- 
ural selection we get, on average, the number of adults with survival probability 
J 2 which is given by 

N 2 (J 2 ,t+l) = J 2 N 1 (J 1 ,t). (4) 

Each one of these adults will generate m 2 descendants with survival probability 
Jo since the mutations are not inherited. 

In general, the number of babies with survival probability Jo is given by 

N Q {J ,t) = m 1 N l {J u t)+m 2 N 2 {J 2 ,t), fort>l (5) 

where J\ — aJo and J 2 — 0(1 — Jo) 1//x . If we substitute equation (5) into (3) 
we can write the following recursive matricial equation 

/ N^Jut + l) 
\ N 2 (J 2 ,t+l) 

where A is the matrix 

. _ ( mxJi m 2 J 1 

{ J2 

Iterating the equation above and using the initial condition, we get for t > 2 

(Sj;;:!)^^°)^ 2 (1 Ji )' (7) 



- A { N 2 (J 2 ,t) ) 
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with A meaning the identity matrix. 

The complete dynamics of the Partridge-Barton model can be obtained by 
diagonalizing the matrix A. We have, explicitly (for t > 2) 

Ni(Ji, t) = === mi Ji (A + -A_ ) 

y mj Ji + 4m 2 Ji J 2 

+m 2 J!J 2 (A^ 2 - A'r 2 )] , (8) 

» r .n </i-7V (Jo,0) r , t _ 2 , t _ 2 N 

N 2 (J 2 ,t) = ==[m 1 J 1 J 2 (A + -A_ ) 

ymi + 4m 2 Ji J 2 

+m 2 J 1 J 2 2 (X l + 3 - A'r 3 )] , (9) 



where 



TOiJi ± yJm\Jl + 4m 2 J 2 
A± = ^ ^ 10 J 

are the eigenvalues of the matrix A, Ji = aJ an d </ 2 = — Jo) 1 ^- Let us 
point out that the time evolution of the babies distribution No(Jo,t), can be 
calculated using equations (5), (8) and (9). Having the expressions above, we 
can determine the evolution of many other quantities like the total number of 
inviduals at age i Ni(t) = J Q Ni(Ji,t)dJi or their mean survival probabilities 

f 1 JiNi(Ji,t)dJi 

(Ji)(t) — ji N ^j ^ dJ ■ The given input parameters arc the initial population 
No, the birth rates (mi and m 2 ), the mutation strengths (a and j3) and the 
pleiotropic constraint (x). 

4 Asymptotic Limit 

Before taking the asymptotic limit, we observe that A+ is the largest eigen- 
value for all possible values of the input parameters. Once these parameters are 
fixed and J 2 = (3(1 — Ji/a) 1 ^, A + is in the last instance a function of J\. From 
the equation (8) we have, asymptotically 

JVi(Ji,t) w e *MA+(Ji)]. (11) 

By integrating in Ji the expression above, we can get the total number of 
juveniles N\{t). It is convenient to change the integration variable J\ for a new 
variable y (a monotonically increasing function of Ji), y = —cot(irJi), such that 
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r°° ptln[\+{y)] 

For large time t, this integral can be evaluated by the saddle point approxi- 
mation. We thus obtain 

t ln[X+(y)] 

Ni(t) = A(y) ^ , (13) 



where y is the value which maximize the eigenvalue A+ and A(y) = / ^ 

In the original paper of Partridge and Barton, the optimization process was 
achieved by a direct (and not well explained) maximization of the growth rate. 
Here, in our formalism, it is a simple and a natural consequence of taking the 
asymptotic time limit in the exact evolving equations. Further, the growth rate 
or the Malthusian exponent is simply given by ln[X + (y)). 

To have deepest insight in the dynamics, let us determine the probability 
density Pi ( Ji , t) of finding a juvenile at time t with survival probability between 
Ji and Ji + dJi . It is given by 

P^t) = ^) = ^4 ^Vie tln[ ^ ] (14) 
[ ' Jo Ni(Ji,t)dJi 1 ' 

where we have used equations (11) and (13) and y = —cot(jrJi). Clearly, at 
the asymptotic limit, the distribution probability P\{J\ 1 1 — > oo) approaches the 
Dirac delta function 5(Ji — Ji) and the mean survival probability at age 1, is 
simply given by (Ji) = J\. Similar results can be obtained for the ages and 
2. Another interesting quantity which can be calculated is the population mean 

Y^ 2 i Ni(t) 

age (A) (t) defined as (A) (t) = ^=|f . It is straightforward to show that 

( A )(t) = 7 + 2 + { 27(l + ^i)+ 7(1 + m 2 ) ] ,-1 0( ,- 2) 

WU 7 (l + mi) + (l + m 2 ) \2[ 7 (l + m 1 ) + (l + ™ 2 )] 2 J + 1 ' 

(15) 

where 7 = x+ i Jl > with J 2 = /3(1 — Ji/a) 1 ^. So we rederive, in a quite simple 
way, the power law decayment first found by Ray M]. 
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5 Discussion 



We solved exactly in this paper the Partridge-Barton model under the action 
of arbitrary pleiotropic constraints and deleterious somatic mutations. Through 
a matricial formalism we were able to predict the complete time evolution of 
the population. We derived analytic expressions for the time dependence of the 
mean survival probabilities and the Malthusian exponent. Since for large time 
t the system behavior is controlled by the largest eigenvalue, it was possible to 
obtain the steady state values of the survival probabilities and to demonstrate, 
in a simple way, that the population mean age has a power law i _1 decayment 
to its final constant value. 

For comparison with our analytical results, we also performed some Monte 
Carlo simulations. In these simulations, the natural selection is implemented by 
discarding any individual with survival probability smaller than a random num- 
ber (generated from a uniform distribution). The deleterious somatic mutations 
and the antagonic pleiotropy can be easily incorporated into the computer pro- 
gram. More difficult is to avoid an explosion of the computer's memory due to 
the unlimited growth of the population. To take this problem into account, we 
resort to the Verhulst factor p^ j which is commonly used in such circumstances. 

In Figure 1 we put together the analytical solution and the Monte Carlo re- 
sult. The exact solution was plotted by inserting equations (8, 9 and 10) into the 
expressions for the mean survival probabilities (Jj)(£) and by integrating them 
using the software Maple jlj. We conclude that the Monte Carlo simulations 
confirm very well the theoretical results. 

Finally, let us to point out that, unfortunately, the technique developed here 
cannot be applied to the case in which mutations are hereditary. The main 
reason for this come from the fact that the equation (5) is no longer valid. 
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FIGURE CAPTION 



Figure 1. The continuous lines correspond to the analytical solutions and 
the square symbols to the Monte Carlo simulations. We used a = 0.82, [3 = 0.67 
, x = 4, mi = mi = 1 and No = 4000. The steady state values are J\ = 0.77 
and Ji = 0.33. There is senescence, i. e., J2 < J\. 
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